Phase transitions for a collective coordinate coupled to Luttinger liquids 
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We study various realizations of collective coordinates, e.g. the position of a particle, the charge 
of a Coulomb box or the phase of a Bose or a superconducting condensate, coupled to Luttinger 
liquids (LL) with N flavors. We find that for Luttinger parameter | < K < 1 there is a phase 
transition from a delocalized phase into a phase with a periodic potential at strong coupling. In the 
delocalized phase the dynamics is dominated by an effective mass, i.e. diffusive in imaginary time, 
while on the transition line it becomes dissipative. At K = | there is an additional transition into 
a localized phase with no diffusion at zero temperature. 



Diffusion and propagation of massive particles sur- 
rounded by a bath is one of very challenging problem 
of condensed matter. Historically it started with the cel- 
ebrated Brownian motion pQ in which the interactions 
between a classical particle and the microscopic motion 
of the classical bath, lead to a diffusion, connected by the 
Einstein relation to a finite friction. 

This problem gets incredibly more complicated when 
the bath becomes quantum. Indeed the excitations of 
the bath can lead, by Anderson orthogonality effects to 
a modification of the motion of the quantum particle or 
the collective coordinate coupled to the bath [2]. One of 
the realization of such problem is the polaron problem 
[3] where the interaction with the vibrations of the lat- 
tice leads to an increase of the mass of the particle and 
even potentially to self trapping. This type of problem 
has recently benefitted from the recent progress in cold 
atomic systems [I]. Indeed in such systems impurities in 
quantum baths can be realized in a variety of manners 
ranging from fcrmi- or bose- mixtures to ions in conden- 
sates, and at various dimensionalities [5Tll4|. 

A situation of special interest is provided by a one- 
dimensional bath for which the bath-bath correlations 
can become highly non-universal due to the existence of 
the powerlaw correlations characteristics of a Luttinger 
liquid (LL) [TS]. In that case special effects can poten- 
tially occur, as clear from the static impurity case [16] 
and mobile ones coupled to single baths [T7] [TS]. In 
particular it was shown recently [19] that this led to a 
new universality class for the motion of the impurity, for 
which in particular subdiffusion can occur. This very 
rich situation was explored further. On the theory side 
both diffusive [20H23], kicked [13 [H] and driven impu- 
rities [2"5"H2"7] were considered. On the experimental side 
driven impurities [12], mixtures of 87 Rb and 41 K [13] and 
87 Rb experiments with local addressability [14] were suc- 
cessful implementation of the one dimensional problem. 



In this paper we study the physics of a collective co- 
ordinate coupled to N Luttinger baths (N ^> 1), e.g. a 
particle position, charge of a Coulomb box, a phase of a 
Bose-Einstein condensate (BEC), or a phase of a super- 
conducting grain. For concreteness, the presentation uses 
the particle coordinate language while other realizations 
are discussed below. We show that this problem exhibits 
a novel localization-delocalization transition as a func- 
tion of the bath characteristics, as summarized in Fig. [T] 
The localized and delocalized states are separated by a 
line on which the motion is simply diffusive. We discuss 
the consequences for the linear response of the impurity 
to an external force, and propose potential experimen- 
tal realizations of this device in cold atomic systems and 
condensed matter ones. 

We consider a particle of mass M coupled to a LL with 
a contact interaction H mt — gp(X) where X is the opera- 
tor measuring the impurity position and p(x) the density 
is the LL. We study this model in the large N limit, 
so the impurity becomes coupled to N independent LL 
and the interaction becomes Hi nt = g X^^Li Pi G^Q ■ The 
action of the system can be computed by a cumulant 
expansion in powers of g and only the second order cu- 
mulant remains when g 2 N = 0(1). Indeed the fourth 
order cumulant is of order g 4 N ~ 1/N and can be ne- 
glected. Using the expression of the density in a LL [15] 
P (x,t) = P u-\d x (j){x 1 T)+pv[e 2l{7TpoX " t ' {x > T ' )) +h.c.] where 
4>{x, t) is the bosonic phase, and performing the Gaussian 
integration over the LL Hamiltonian, the action becomes 
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where we have used the dimensionless variables X = 
2irp Q X and M = M/(2irp„) 2 , r) = 2tt g 2 p^N '/ 'A 2 : , r is the 
imaginary time, u the velocity of excitation in the LL, 
and K the Luttinger parameter that controls the power- 
law decay of the correlation functions. A frequency cutoff 
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A = u/ct is used to have a dimensionless coupling 77 where 
a ~ 1/po is the natural momentum cutoff of the LL. In 
the above expression only the oscillating (backscatter- 
ing) term in the density has been retained. Indeed the 
d x <p(x,r) interaction can be integrated, leading at long 
times to (X T — X t i) 2 /(t — r') 4 , i.e. an cu 3, term in fre- 
quency which can be neglected relative to the bare ki- 
netic energy term of the impurity Muj 2 . We have used 
that for a LL one obtains [13] (^4>{X T ,r)-i^X T ,,r')^ LL _ 

[(X(t) - X{t')) 2 + u 2 {t - t') 2 ]- k . We have also made 
the additional assumption, which will be verified in what 
follows that the impurity is less than ballistic and thus 
that ((X T - X T ,) 2 ) < u 2 (t - t') 2 . 

To solve for the thermodynamics of we consider 
first a renormalization group (RG) process [55] valid for 
large rj, which was also applied to the K = 1 case [2"D] . 
The action ([I]) is approximated by its short time form 
where it becomes Gaussian 
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where J T 1 ~ c t ° 2 T t) = -2L(1 - 2K) sin(Jfrr) 
hC k \lo\ 2K - iT so that C K = l- 0.85(K - 1) + 0(K - l) 2 . 
The cutoff A is replaced by A' and the interaction is aver- 
aged with So in the small frequency interval A' < to < A 
leading to dA/irCKfjA. The action has then a renormal- 
ized coefficient rj R (A') 2 ~ 2K where 
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with Ck — > 1 to 1st order in either I/77 or 1 — K. Hence 
if K > 1 j] R flows to small values, while if K < 1 there 
is an unstable fixed point at r\ c = 2 tt(i-k) ■ V > Vc flows 
to large values, while 77 < r\ c flows to smaller values of 
rj. One can integrate ([3]) when i] < rj c down to rj R ~ 1 
below which the RG is not controlled. The new cutoff is 
interpreted as an effective mass [29, 30 M* 
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which for 7rr;(2 — 2K) <C 1 but nr) ^S> 1, i.e. far from the 
transition point, represents an exponentially large mass 
M* ~ e 7 ") as for the K = 1 case [13 [3D]. 

To confirm this scenario, and study the properties of 
the three resulting phases, we follow a variational scheme 
[31j where we find the best quadratic action approximat- 
ing the original action ([I]). The corresponding Green's 
function l//(w) is a solution of the self-consistent equa- 
tion 
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We note first that at ui w A the solution is /(u)-Mw 2 
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In the following we focus on w <C A and 



on if < 1. As a first option we consider f(uj) — 
ri* Cku 2K ~ x / k 2K ~ 2 . The integral in the exponent con- 
verges as r — > 00, so it is J Q A dw'/ (7rf(u)')) = [irrj* Ck(2 — 
2K)]~ 1 , hence (pM) reduces to 
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This equation has solutions only if 77 is sufficiently large, 
i.e. ttCk ??(2 — 2if ) > e. A second possible solution is 
f(u) = Tj*\oj\. The exponent behaves as J Q A 1 ~ C °, S 1 ^ T = 
^rlnAr, since the r integral is dominated by long r, 
hence 
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which is a consistent solution on a line — = 2(1 — 
If). The third possible solution is similar to the bare 
one f(co) = M*w 2 , then the exponent behaves as 
Jo" 1 ™? = M/2M*, leading to 
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where the r integral is cutoff by M* and then an expan- 
sion in ui is possible. This solution is acceptable only for 
intermediate or weak coupling and then M* ps M. 

The existence of the three possible regimes in the varia- 
tional approach shows that there is a connection between 
this method and the RG process [3T]. One can replace 
(§by 
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where B(rj) can be found as a perturbative series. RG 
proceeds by defining scaling functions such that the 
frequency dependence is replaced by a renormalized 
fj(id), i.e. f(id) = fj(id)id with 77(A) = 77, f'(u) = 
B[fj{u))}, f»(u) = £C[fj(u)]. The flow is given by 



dfi(u) , C(rj) 

| A = B{ri) - 77 = 
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i.e. B{ri) = rj and C(rj) = are fixed point conditions. 
Eq. ([9]) and its derivative at id — A yields 



B'{rj)[B{n)-r l ][-B l {n) + 2K- 



irrj 



(2K - 2 + — )B{rf) + B"(rj)[B(rj) - r/} 2 (11) 

■KTj 

with boundary conditions at a large 770 from a perturba- 
tive analysis B(r) ) = Tj (2K - 1) + i, £'(770) = 2if - 1. 

Note first that a fixed point .8(77) = 77 is achieved at 
2K — 2 + — = to all orders, consistent with Eq. 

7TTI 5 ^ 

as an unstable fixed line. Note next that f{td) ~ kT is 
a solution with B(rj) — 2?j, consistent with Eq. It 
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corresponds to a fixed point at fj(ui) ~ cj —> 0. Finally, 
for large 77 one can solve Eq. ^ perturbatively, leading 
to f(ui) — r/ R (uj)uj with j] R (oj) given formally by Eq. (13j) 
with the replacement A' — > uj. Using the RG structure 
this can be integrated and for "C 1 — K it yields 
f(u) = ^A 2 " 2 *^ 2 *^ 1 , i.e. So of Eq. (|} is a fixed point 
action for large 77. 

The analysis leads thus to three different possible be- 
haviors for the system: 

i) At 1 — K = 5^ the system has a dissipative behavior 
with f(cj) = T]\lj\. 

ii) The case 1 — K < 7^— flows to small ri and eventu- 
ally to an M*u 2 form, with ((X T ~ X ) 2 ) ~ |t|, which 
corresponds to a delocalized phase. The effective mass 
M* is identified by the RG flow, as in Eq. Q. Note 
that even in this delocalized phase, some effects of the 
underlying quasi-long range periodicity of the LL with 
the wavevector 2wpo are still felt by the particle. Indeed 
its correlation at that periodicity are only very slowly de- 
caying (cosA r cosA ) ~ t~ 2K . This indicates that the 
particle has a much greater chance to be found at some 
particular places on the chain. This can be understood 
qualitatively by the argument that the particle moves in 
the "charge density wave" of wavevector 2irpo provided 
by the LL, hence the particle diffuses predominantly by 
tunneling between lattice sites spaced by l/po- On the 
mathematical side, this property which is apparent in a 
first order calculation in 77 |31| is in fact known in general 
in the context of XY models with long-range interactions 



iii) The case 1 — K > ^L- flows to large 77 with eventually 
/(w) ~ lu 2K ^ 1 . From this form one could naively expect 
that the correlations of ([X T — Xq] 2 ) to be convergent 
and thus this phase to be a localized one. The situation 
is in fact more subtle and we discuss this phase in more 
details below. 

A summary of the various regimes can be found in 
Fig. [T] and the corresponding correlation functions are 
indicated in Table fl] 

Let us now discuss in more details the physics of the 
phases, in order to ascertain that they are not mere 
artefact of the variational approach. A priori the vari- 
ational scheme is reliable for large 77, i.e. for the tran- 
sition line near K = 1. For K < 1 one can com- 
plement the above analysis by a mean-field approach 
similar to the one used in the context of XY mod- 
els with long-range interactions [33] ■ We take h = 
(cos X T ) as an order parameter. The interaction term 
in Q decouples as r]A 2 ~ 2K h J t cosX t J t , \t - t'\~ 2K = 
r]Ah 2I ^_ 1 J r cosX T . The self consistency relation, linear 
in h, is 1 = r/A^ry f T , (cos r cos T /) = 4j)MA^j; it 
yields the critical line rj c — 2 tttx above which (cos X T ) ^ 
0. We expect the mean field result to be more reliable 
near K = |, where the range of the interaction increases. 
As K increases from K = | fluctuations will increase the 
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FIG. 1. Phase diagram for an impurity in a bath of LL 
as a function of the LL parameter K and the interaction pa- 
rameter between the impurity and the bath r\. Four regimes 
can occur (see text) in which the impurity is delocalized, just 
dissipative, periodically localized, or localized. Dashed lines 
indicates boundaries out of the control of the perturbative 
RG. The corresponding correlation functions are given in Ta- 
bled 



critical value, eventually joining the transition line with 
the variational form r/ c — 27r (i-K) near K = X. Note that 
mean field exponents become valid [33] when K < 3/4, 
e.g. (cos X T ) ~ ^/r\ — t] c . In Fig. [I] we plot the transi- 
tion line as an interpolation between the mean field at 
K < 3/4 and the variational form at K > 3/4. We see 
that the point K = 1/2 plays an important role, not 
captured by the variational or RG approaches. Below 
this point the interaction is so-long range that an or- 
dered phase would exist, within the mean-field solution, 
for arbitrary strength of the coupling 77. 

We reconsider now the periodic phase. Since the or- 
dered phase would have many degenerate minima of the 
order parameter, we expect the presence of instantons. 
Such instantons are known for the K = 1 case [541, [35] . 
To estimate the action of such instantons we assume an 
infinitesimal ordering field that selects the ground state 
and consider then interaction energies as in (JlJ. Assum- 
ing an instanton with width To, the interaction term has 
the form t](At ) 2 ^ 2K while the mass term is ~ M/t , 
hence the action is minimized at At ~ ( (i-KyqB K ) 3 ~ 2K 
for K < 1; the numerical prefactor Bk is known at 
K = 1, Bi = 7r. Note that the mass term and K 7^ 1 set 
a finite scale for tq , unlike the K = 1 case. The instanton 
action is then 



MA 



(2 - 2K) V B K 
MA 



3-2K 
2-2K 



(12) 



Such instantons mean that the coordinate X T can tun- 
nel between neighboring minima of the ordered (cos X T ). 
Assuming independent instantons this would imply that 
((X T — Xq) 2 ) = D\t\ has a finite diffusion constant. The 



4 



TABLE I. Correlations of the phases in Fig. 1 at T = 0. 



correlation 


delocalized 


dissipative 


periodic 


localized 


(COS Xr) 








constant 


1 


(cos X T cos Xo) 


~ \r\~ 2K 


^ | r |-(2-2K) 


constant 


1 


((X T - Xo) 2 ) 


~ r 


~ In t 








tunneling rate and D are then proportional to the instan- 
ton density, i.e. D ~ e -Si " st . 

In particular we consider K — > \ and an instanton 
localized at r = 0. The dominant contribution for the 
instanton center at |r| < t comes from |t'| > r that 

which 



involves \X T \ > \X T , \ and jj T ,| >To 
diverges at K 



-2K 



2K-1 



|, hence 



\M I (Xf T 



7]A(At ) 



tt(2K 

S' comes from the instanton tails where X T , X T < are small 

(up to 2tt) and comparable. This action is similar to 

the well known sine-Gordon system, identifying Bk ~ 

(2K — whose instanton (or soliton) solution has a 

i _ i 

width t ~ (2K — l)a and action S inst ~ (2if - 1) 2. 

Assuming independent instantons the diffusion constant 
would diverges at X = §, i.e. lnZ3 ~ (2 A" — 1)~2. Wc 
propose that the whole range of the periodic phase in Fig. 
1 has instanton solutions with a finite action, with an 
explicit solution provided by the sine-Gordon system at 
K — > 1/2. The tail contribution S' has a structure close 
to that of the 2nd term of Eq. (J2j|, that may affect the 
interaction between instantons. Therefore, given the long 
range form of the interaction, to ascertain the correct 
behavior of ((X T — Xo) 2 ) at large time requires further 
study of how these instantons interact, which is left for 
the future. 

We consider next the system at K < 1/2. This case 
has been studied in the context of discrete XY models 
[36H38] and was shown to have a phase transition in the 
limit that the coupling vanishes as a power of the system 
size, which in our case is /3 = 1/T, i.e. there is a crit- 
ical value for rj(f3) 1 ~ 2K . Hence at T = the system is 
fully ordered and cosX T = 1. Furthermore, instanton ex- 
citations would involve the effective coupling ri(f3) 1 ~ 2K , 
hence will have diverging action. Below a putative mean 
field transition temperature T x where Ar/MAT^^ 1 = 1 
the system dynamics would be similar to that of the pe- 
riodic phase. 

Let us conclude this part by noting that the hypothesis 
made at the beginning to neglect X T — Xo compared to 
r is indeed justified in all the phases. Furthermore, note 
that although the results of the present paper are derived 
in the large- X limit, we of course expect them to extend 



to a finite number of component as well. E.g. for the 
Coulomb box case deviations due to finite N appear at 
exponentially small temperatures |39j . 

Finally we discuss possible realizations of our model 
with various collective coordinates that are potential can- 
didates for experimental studies. A first possibility to ob- 
tain ([lj is to use a fermion Coulomb box [3U]. Following 
the Ambegaokar-Eckern-Schon mapping |41j one intro- 
duces a phase X T such that X r measures the charge in 
the box while the charging energy corresponds to 1/M. 



f(l -cosX T ) + S' (13) Thekernelin Q is then £ Q G q ,^t - £ fe G M (r' - 
~ 1) J T where i is the channel index, a, k are internal qi 



where i is the channel index, a, k are internal quan- 
tum numbers of the dot and LL, respectively, and the 
Green's functions are for either free fermions on the 
dot, ~ 1/(t — t') while for fermions in the LL it is 

~ |r — T '\~2( Kf+1 / Kf K Hence an effective action of the 
form (TT1) with 2K = 1 + \{Kj + l/Kf), realizing only 
K > 1 cases. In case that the LL terminate at the 
Coulomb box a boundary Green's function [IS] is needed 
G x= o,i(t - r') ~ \t - r'|- 1/K/ , hence 2I\ 1 • 
and K < 1 is realizable. A second possible realiza- 
tion corresponds to a BEC with a phase 9 T that weakly 
couples to bosonic LL's with boson operators ^(t) as 
ge Tl I r „(T) + h.c The average involves now the boson's 
Greens function ~ |r — r'| 1 / 2/fb , — > 00 for noninter- 
acting bosons and Kb decreases to 1 for on site repul- 
sion U — > 00. Hence ([T]) is realized with K = 1/4 A), 
and K < 1 is possible. In analogy with BEC, a super- 
conducting grain can Josephson couple to fermionic LL's 
providing a third type of realization; here again K < 1 
can be realized. 

In conclusion, we have studied the physics of LL envi- 
ronments that couple to a collective coordinate such as a 
impurity position, charge of a Coulomb box, a phase of a 
BEC or that of a superconducting grain. We have shown 
that the coupling to the bath leads to various phases for 
the collective coordinate ranging from delocalized, dis- 
sipative, periodic and localized. Our results are sum- 
marized in Fig. [T] and Table |TJ showing the distinctions 
among the various phases. We believe that the large set 
of realizations for the collective coordinate and the vari- 
ous phase transitions will stimulate further research. 
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